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We perform a detailed study of a specific Two Higgs Doublet Model (2HDM) with a U(l) gauge 
symmetry, instead of a typical Z2 discrete symmetry, containing a very light gauge boson Z' (GeV 
scale or below). The Standard Model (SM) fermions do not carry U(l) charges, but induced cou- 
plings to the Z' (called the dark Z) are generated through mixing with the SM neutral gauge bosons. 
Such a light Z' could explain some astrophysical anomalies as well as the muon g — 2 deviation, and 
has been the subject of great experimental interest. We consider the scenario in which the 125 GeV 
SM-like Higgs (H) is the heavier scalar state, and focus on the lighter neutral state (h) as well 
CO ■ as charged Higgs. We analyze the constraints on the model from various experiments and predict 

novel channels to search for these Higgs scalars at the LHC. In particular, experiments looking for 
, lepton-jets are among potentially important searches. 

CN 

Ph. I. INTRODUCTION 

< 

. At long last, the experimental exploration of the Higgs sector of the electroweak theory has begun. The discovery 
of a scalar boson with a mass of approximately 125 GeV at the Large Hadron Collider (LHC) is a spectacular triumph 
of the Standard Model (SM) [H,[|. The hints of a larger than expected branching fraction for H — > 77, although they 
have recently weakened, could suggest possible extensions of the scalar sector Q or of other sectors. (For some vector- 
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like fermion extension examples, see Refs. [I|-Q1|-) Certainly, detailed analyses of these extensions are worthwhile and 
C^H' might provide guidance to experimenters searching for physics beyond the SM. 

The simplest and most studied such extension is the Two Higgs Doublet Model (2HDM), recently reviewed in 
t" 1 h Ref. A potential problem for such models are tree-level flavor changing neutral currents (FCNC) which will 

occur whenever fermions of a given charge couple to more than one Higgs multiple t fl^. [l8j . While it is possible 
to ameliorate the constraints by assumptions about the flavor-changing couplings |l9l. |2Q||. the usual solution to 
the problem is eliminating the tree-level FCNC altogether by using a discrete Z2 symmetry. Depending on the 
transformation properties of the right-handed fermions under the Z2, several different models can be obtained [l6| . 
The most familiar two models are the Type I model, in which all of the fermions couple to a single Higgs multiplet, 
\Q \ and the Type II model, in which the down-type quarks and charged leptons couple to one Higgs multiplet and the 
up-type quarks couple to the other. 

A Z2 symmetry can be promoted to a gauged U(l) symmetry as discussed in Refs. [2l], [22|. In Ref. [2l|, the SM 
fermions do not have charges of the new U(l) gauge symmetry and a Type I 2HDM model results. (See Ref. [23[ 
for another example in the fourth generation context.) In Ref. [22j], the SM fermions have charges under the U(l) 
gauge symmetry, resulting in various types of 2HDMs. As emphasized in Ref. 12211 . one of the interesting features in 
this kind of model is the absence of the pseudoscalar A. It is a common effect [24| of a spontaneously broken U(l) 
gauge symmetry that the pseudoscalar is eaten to become the longitudinal component of the new gauge boson Z' ' . Of 
course, one could restore the pseudoscalar by adding a complex Higgs singlet to the model. 

In our paper, we will closely follow the scenario of Ref. [21]. In this model, all SM fermions are neutral under the 
new Abelian gauge group U(l)', yet the Z' coupling to the SM fermions can be mediated by mixing with the SM 
neutral gauge bosons. The Z' can communicate with the SM fermions only through kinetic and/or mass mixing with 
the SM neutral gauge bosons. The fact that the Z' does not directly interact with fermions eliminates the numerous 
constraints on Z' masses, and allows the Z' to be extremely light, with mass allowed below 0(1) GeV. (For a review 
on a complementary heavy Z' , see Ref. [25j-) 

This model in Ref. [2l| was not originally proposed in the context of 2HDMs. It was first proposed as a generalization 
of the so-called dark photon model whose coupling is of the same form as the photon coupling. While the two models 
both adopt kinetic mixing of the U(1)y and the U{1)' , the difference comes from the fact that the Z' gets its mass 
only from a Higgs singlet in the dark photon model while it gets the mass from a Higgs doublet (and also from a Higgs 
singlet if a Higgs singlet exists) in the dark Z model. As a result, the dark Z can couple to the weak neutral current 
as well. When the Z' is very light, it can explain [26-28] the astrophysical anomalies such as the 511 keV gamma-ray 
from the galactic center observed by the INTEGRAL satellite 29] or the positron excess observed by ATIC and 
PAMELA experiments [3l[, depending on the property of the dark matter. It can also explain the 3.6a deviation 
of the measured muong — 2 from the SM prediction [32|, |33[ . (For some discussions about non- Abelian dark gauge 
sectors, see Refs. [33-l37|. See also Ref. [38[ for a discussion on W-W mixing.) 
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It is remarkable that a very light gauge boson, introduced to provide a dark matter explanation to some astrophysical 
anomalies, can also provide an explanation for the absence of tree-level FCNC in a 2HDM. It leads to a very different 
phenomenology from the ordinary 2HDMs. Some of the implications for the dark Z boson as well as the 125 GeV 
SM-like Higgs state, H, of this model have been discussed in Refs. [U H^, H(| ■ In this paper, we will mainly discuss 
the phenomenology of the other Higgs bosons of the model, and refer to this model as the "Dark 2HDM" to emphasize 
that our study of this model focuses on the comparison to the ordinary 2HDMs. We particularly focus on the scenario 



that the other Higgs, h, is lighter than the SM-like Higgs, H, of 125 GeV, which was a subject of the study in Ref. [41 1 
for the ordinary 2HDMs. We will see that various constraints force one to add a singlet to the model. The decay of 
the h to Z'Z' will provide a dramatic signature that could be on the verge of discovery at the LHC, and the decay of 
the charged Higgs to hW —> Z'Z'W would also lead to striking signatures. 

In Sec. |TT1 we describe the model and the particle masses and couplings. In Sec. IIII1 we consider the constraints 
from LEP and the LHC on the neutral Higgs bosons, and in Sec. IIVI study the remarkable signature of the light Higgs 
and the implications for LHC experiments looking for multi-lepton jets. In Sec. [V] we discuss the properties of the 
charged Higgs boson and in Sec. IVI| we present our conclusions. A discussion of the properties of the Z' boson is in 
Appendix. 



II. THE MODEL 

In this section, we describe the Dark 2HDM proposed by Davoudiasl, Lee, and Marciano [2l|. The scalar sector is 
composed of two doublets and a singlet under the SU(2)l. As we will see, the model with just two doublets and no 
singlet is not compatible with various constraints. Including a singlet can avoid this issue. 

The gauge group is SU(3)c x SU(2)l x U(1)y x U(l)', and the SM fermions are assumed to be neutral under the 
extra U{1)'. The gauge part of the kinetic Lagrangian includes kinetic mixing between the U(1)y and U(l)' [42[ and 
is given by 

£ ~ = -\ B ^ y + \^e^ B ^ z '" u \V" V (!) 

where B^ v = d^By — d v B^ and Z' = d^Z^ — d^Z'^. The hat notation indicates the states before diagonalization. As 
discussed in various places, for examples see Refs. [U H^|, one can redefine the fields to remove the kinetic mixing 
term parametrized by e, which is experimentally constrained to be small. 

When there is one Higgs doublet and the Z' gets the mass from a Higgs singlet, this leads to an induced coupling 
of the Z' to the electromagnetic current 

Aiark 7 = -eeJ^ {J? m = Qff-ff) (2) 

and the induced coupling of the Z' to the weak neutral current is negligibly small because of the suppression by high 
orders of e. This gauge boson is widely called the "dark photon", since fermions couple to the Z' with a very small 
effective electric charge (eQf) for a given electric charge (Qf) of the fermion /. 

The model in Ref. [2l[ added an additional Higgs doublet charged by the U(l)' which introduces mixing in the 
Z-Z' mass matrix. The Z-Z' mass mixing, parametrized by ez, is also taken to be small. The value of ez will be 
determined by parameters in the scalar potential. Because of this mixing, the Z' couples to both electromagnetic 
current and weak neutral current as 

£daxk z = - (eeJ» m + ezgzJ^c) % (j£ c = [±T 3/ - Qf sin 2 6 w ]f^f - [^T 3/ ]/V 7s A (3) 

with gz = gj cosQw — g'/ s'm9w ~ 0.74, T 3 f = ±1/2 and the weak mixing angle sin 2 9w ~ 0.23 . The Z' with more 
general coupling (compared to the dark photon) in this model was named the "dark Z" to emphasize that it couples 
to the weak neutral current. In particular, when the e = limit is taken, the Z' couples only to the weak neutral 
current as the SM Z boson does with a suppressed coupling (szgz)- The definition of Ez and constraints on m^, e, 
ez can be found in Appendix. 

We use a similar but slightly different notation from Ref. [111. 1 In conventional 2HDMs, one defines $2 to be the 
doublet that couples to the top quark (or, in the Type I 2HDM, to all fermions). Following this notation, we assume 



1 In Ref. [2ll |. the field $1 was chosen to couple to the SM fermions rather than <J>2, with definitions of tan/3 = v%/vi, tan/3d = V2/vs- 
This leads to different typical values of tan and tan fi^. Care must be taken in defining the heavy Higgs vs. the light Higgs, which also 
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the 17(1)' charges Q'[$i] = Q'[$ s ] = 1, Q'[$a] = with the hypercharges Y[$i] = Y[$ 2 ] = 1/2, Y[$ s ] = 0. The 
is spontaneously broken when <£>! or $5 gets a vacuum expectation value (vev). 
The scalar part of the kinetic Lagrangian is given by 

Scalar = | | 2 + | £> M $2 | 2 + | $ S | 2 (4) 

— ^-mass ~T~ -^coupling H~ * * * (5) 



where 



= (d,, + ig'Y^B,, + igT 3 [<Z>i]W 3lJt + igz'tfl^zf) (6) 



and two doublets $i = ( i , \ ) with vevs Vi (i — 1,2), and a singlet $s = -jr;{vs + 0s + * r /s) with a 

vev vs- The U5 = limit corresponds to the doublets only case. 

The Z and Z' masses, with v 2 — v\ + v\ ~ (246 GeV) 2 and tan/3 = v%jv\ y tan (3d = vs/vi are given by 



where 



C mass = \m 2 zo Z°Z a - tfZ°Z'° + ±m z , Z' Z' + ■■■ (7) 



2 1 2 2 / \ 

m zo = ^9z v , (8) 
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|,„ = g% t ( v 2 C o S 2p + vl) + ^— g z ,g>v 2 cot? P + l (-^-) ff'V, (9) 

COS 0w 4 \ COS / 

A 2 = \gz>gzv 2 cos 2 (3 + - ^— gzg'v 2 . (10) 

Z 4 COS (714/ 



This leads to the Z-Z' mixing angle (£) as 



2A 2 

tan2 £=^2 =3- ( U ) 



Z \ _( cos£ -sin£ \ / Z° 
Z' I ~ [ sin£ cos£, I { Z'° 



(12) 



where Z and Z' are the mass eigenstates. The Z-Z' mixing angle is constrained to be very small (|£| <; few xlO -3 ) 
by precision Z pole measurement at LEP [33| . 

We are primarily interested in a very light Z' [mz> <J 0(1) GeV]. We will work in the m 2 z , -C m 2 z0 limit throughout 
this paper, which requires g 2 Z iiv\ + Vg) <C \g z v 2 as |e| is very small. In this limit we have 
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(13) 

3 2 ,(z, 2 cos 2 sin 2 P + v%) = g 2 ,v 2 ^JL(l - cos 2 f3 cos 2 (14) 

cos" 1 /id 



(15) 

The Z' approaches the massless limit as g Z i — > or vi, Vs — > 0. With the <5 notation defined in Appendix , we can 

S 2 1 
cos a 3 



use cos 2 fid ~ — W (the doublets only limit corresponds to <5tan/3 ~ 1), and have 



g Z 'V cos 2 j3 

mz> ^ j ■ (16) 



amounts to a definition of a. We provide a comparison of notations in the following. Ref. [2l| : <t>i couples to the SM fermions, which 
requires typically small tan/3 (as vi ~ u). o = corresponds to the heavy Higgs ~ SM-like Higgs limit. This paper: <t>2 couples to the 
SM fermions, which requires typically large tan/3 (as V2 ~ i 1 ). a = ±7r/2 corresponds to the heavy Higgs RJ SM-like Higgs limit, (a = 
corresponds to the light Higgs SM-like Higgs limit.) The Higgs to gauge boson couplings in Eqs. I|18[ l - (1 23 1) can be read for Ref. [2~D | 
by cose* +-> sin a, cos /3 -o sin/3 (or iq f> U2), and overal sign flips for C^zz, C^zz'y Chz' z' ■ 
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For simplicity, we assume no mixing between the doublets and singlet and allow mixing only between the doublets. 
The pure doublets case can be reached if we take v,s = (corresponding to cos f3 d = 1) in Eqs. © and (JH) , an d it 
would not change the following couplings in Eqs. (|18|) - (f23| as well as the approximations in Eqs. (|24|) - (f29)) . 

The relevant couplings of vector bosons to heavy Higgs (H) and light Higgs (h), with no mixing between the 
doublets and singlet, are 

•^coupling — -^ChzzHZZ + Chzz'HZZ' + -Chz' z' H Z' Z' 

+ ]^C h zzhZZ + C hZZ >hZZ' + hz hZ . z >hZ'Z' + ■■■ (17) 

where 

Chzz = Crzz (cos /3 cos a [cos£ + (2gz>/gz + £ tan#w) sin£] 2 + sin (3 sin a [cos£ + e taiidyj sin£] 2 ^ (18) 

Chzz' = C H ™ Z (cos /3 cos a [cos £ + (2gz>/gz + etan0 w ) sin£] [sin£ - (2g z >/gz + £ taxi6 w ) cos£] 

+ sin /3 sin a [cos £ + s tan 8\y sin £] [sin £ — s tan #w cos £]) (19) 

Chz'Z' = Chzz fcos/3cosa [sin£ — (2gz'/gz + £tan6*w) cos£] 2 + sin/3 sina [sin£ — etanfl^ cos£] 2 ^ (20) 

and 

Cfczz = Cfrzz (sin/3 cosa [cos£ + £ taixOw sin£] 2 — cos /3 sin a [cos£ + (2gz>/gz + £ tan#iy) sin£] 2 ^ (21) 

Cftzz' = Cf^ z (sin/3 cosa [cos£ + etanfl^ sin£] [sin£ — etan#w cos£] 

— cos /3 sin a [cos£ + (2g^/ jgz + et&wOw) sm£] [sin£ — (2gz>/gz + £ tanflwO cos£]) (22) 
ChZ'Z' — Cf^zz (sin/3 cosa [sin£ — e tan 6>vy cos£] 2 — cos/? sin a [sin£ — (2gz' jgz + et&nOw) cos£] 2 ^ (23) 

with the SM Higgs-Z-Z coupling = Note that, unlike Ref. [2l[, we include couplings with a general 

mixing angle a. 

Since |£|, |e| <C 1, we can get approximations, using Eq. ([T5"j) . 



C/zzz Cf^cos(/3-a) (24) 
C^zz' - -Cf^ z (2 3z V. 9z )cos/3sin/3sin(/3-a) (25) 
C H z>z> Cf| z (2 5z V. 9z ) 2 cos/3sin/3(cos 3 /3 sin a + sin 3 /3 cos a) (26) 



and 



C hZZ =s Cf| z sin(/3 - a) (27) 
CW< =s Cf| z (2. 9z V 5z )cos/3sin/3cos(/3-a) (28) 
Cfcz'*' - Cf| z (2. 9z V. gz ) 2 cos/3sin/3(cos 3 /3cosa-sin 3 /3sina) (29) 

giving the expected 2HDM couplings of the two neutral Higgs bosons to the Z in Eqs. and (|2"T|) . In the a = tt/2 
limit, we can reproduce the relation shown in Ref. [2l|, i.e. e z ~ Chzz' /Chzz — Chz'Z' /Chzz 1 - 
The scalar potential is 

V = Vi + F 2 (30) 
Vi = m 2 1 $l$ 1 +TO 2 2 ^$ 2 + y($ t 1 $ 1 ) 2 + y(^$ 2 ) 2 + A 3 ($ t 1 $ 1 )(^$ 2 ) + A 4 ($ t 1 $ 2 )(^$ 1 ) (31) 

V 2 = m 2 3 ^$ s + ^($^ s ) 2 (32) 

with all coefficients real. Interestingly, the terms m 2 2 ($1$ 2 + $2*1) and ^[($I*2) 2 + (*2 $ i) 2 ] whose coeffi- 
cients are generally complex are forbidden by the U(l)' gauge symmetry. We assume the terms A7(${<I , i)(<I > ! ? <I , s) + 
As( < I , 2 < J > 2)( < l ) 5 ( i > s) are negligible (no mixing between the doublets and singlet). 
The charged Higgs mass is given by 

2 ^4 2 /oo\ 

m H± = -—v (33) 
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which requires A4 < 0. 

The Higgs singlet mass is given by 

m% = X 6 v 2 s (34) 
and the neutral doublet Higgs mass-squared matrix is given by 

M 2 _ ( Aiuf (A3 + A 4 )uiU 2 \ r orS 

^Higgs - ^ (A 3 + \ A ) Vl V 2 \ 2 v\ ) ■ [6b > 

The mass eigenstates of the doublets are H and h (with m# > rrih) with 

m 2 H = \ [\yv\ + \ 2 v\ + yJJXivf - \ 2 v\f + 4(A 3 + A 4 ) 2 ^^ 
rr4 = 5 (a iW ? + \ 2 v 2 2 - ^(X.vi - \ 2 vlf + 4(A 3 + A 4 ) 2 «m) 



(36) 
(37) 



and the mixing angle is given by 

c 

ta.u 9.n> = - , 

Xivf - \ 2 v. 



2(A 3 + A 4 )wiW2 , . 

tan 2a = — — ~ ^ — o - $8) 



H \ ( cos a sin a \ ( 61 



h J \ — sin a cos a 



>2 



(39) 



Since only $2 couples to the SM fermions in this model, the sin a s» ±1 limit provides the heavier Higgs (H) as the 
SM-like Higgs (H ~ 02 7 ^ ~ 4>i up to a sign), which is the case of interest in the subsequent sections in this paper. 
The other limit, sin a « 0, would have provided the light Higgs (h) as the SM-like Higgs (H ~ <p±, h ~ 02 up to a 
sign). Note that the SM-like limit for the H is achieved for p = a instead of j3 — a + it/2 which would be true when 
the light Higgs h is the SM-like one. 

The relative coupling of the Higgses to the SM fermions and W boson in the Dark 2HDM is the same as the Type 
I 2HDM as following. 

S1H Oi 

Htt, Hbb, Htt : -, HWW : cos(/3 - a) (40) 

sin p 

COS Ot 

htt, hbb, hrr : , hWW : sinCS - a) (41) 

sinp 

The relative coupling of HZZ (hZZ) is the same as that of the HWW (hWW) to a good approximation since the 
Z mixing is very small [Eq. ([2"g)l. Eq. (|2"7)]. 

III. LEP AND LHC CONSTRAINTS ON THE NEUTRAL HIGGS BOSONS, h AND H 

In this section, we consider various constraints/predictions on the neutral Higgs bosons in the Dark 2HDM. They 
include (i) LEP limits on associated production of the light Higgs (Z* — > Zh), (ii) LEP limits on the marginal Z width 
for Z — > hZ' decay, (iii) LHC measured branching ratios in various modes (H -> bb, r+r~, WW, ZZ) including the 
most precisely measured H — > 77. 

We will focus on the case in which the mass of the light Higgs is ran/ 2 <J rrih <J mz, i.e. roughly, rrih — 60 — 90 GeV 
range. If it is lighter than mjy/2, then the fact that H — > hh decays do not dominate restricts sin a to be infinitesimally 
close to sin a = 1 [41]. If it is heavier than mz, then the h — > ZZ' channel opens up although our results may not 
change qualitatively. 

A. LEP bounds from associated production of h 



A SM Higgs of less than 1 14 GeV was ruled out long ago by LEP experiments looking for associated production with 
a Z. In 2HDMs, however, the hWW and hZZ (neglecting small Z mixing) couplings are suppressed by sin(/3 — a), 
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FIG. 1: (a) LEP constraints on the invisible Higgs decay a(Zh)BR(h — > m.v)/a(ZHsM) at 95% CL. It is the same as sin 2 (/3 — a) 
when small Z mixing is neglected, (b) LEP invisibly decaying Higgs bounds (L3) on a and j3 in all 2HDMs including the Dark 
2HDM (when the Z mixing is neglected). The parameter-space within the bands (central part) are excluded, for a given light 
Higgs mass m h = 60 GeV (black), 80 GeV (blue), 100 GeV (red), and 110 GeV (green dashed). The dashed curve (for 110 GeV) 
shows how quickly the LEP bounds converges for mi J> 105 GeV. 



possibly making the h quite gaugephobic. LEP published upper bounds on this factor as a function of the Higgs 
mass, assuming the Higgs would decay into b quarks or r's (4J446J. However, as we will discuss in Sec. IIV1 h can 
mainly decay into Z'Z' in the Dark 2HDM. Since the LEP searches did not cover fermions from a very light Z' 
\mz< <z 0(1) GeV], one must look at LEP bounds on invisible Higgs decays. These bounds at 95% CL from ALEPH 
[47j |. DELPHI [H, L3 0, OPAL HO] collaborations are given in Fig.Q](a). There are also updated ALEPH results, 
which covers the Higgs mass from 95 GeV [Hljj , but this is beyond our range of interest. For simplicity, we will take 
only L3 data, which is roughly close to average values of the LEP data in the nih range we consider in this paper 
(m/j. ~ 60 — 90 GeV). 

For a given h mass, this bound is given as a band in tan/3 — sin a plane in Fig. [TJ (b). The apparent asymmetry 
between the sin a ~ — 1 and sin a ~ 1 regions originates from tan/3 > 0, which makes the region near sin a ~ 1 (or 
a ~ 7r/2) preferable . As the h mass increases from around 105 GeV, the bands quickly converge as the LEP excluded 
region vanishes [Fig.Q](a)]. This is illustrated as a sudden departure of m/j = 110 GeV curve from the = 100 GeV 
curve in Fig. Q] (b) . 



B. LEP bounds from Z -»■ hZ' 



The h can be produced directly from Z decays into an h and a Z', if kinematically allowed. It will appear as a 
contribution to the invisible Z width. 

For a sufficiently light Z', using Eq. (|28|) . we have 

Y(Z -> hZ>) ~ [C hzz ,f-^* (l-^i) (42) 



647tto|, \ m z 

64tt 



(Jtan/3) cos^/3-a) 1 f- . (43) 



As was discussed in Ref. [2lj], the boosted Z' shows the divergent nature as mz> —t due to the enhancement from 
the longitudinal component of the Z' . 

The limit on the new physics contribution to the undetected width of the Z, which can be obtained by subtracting 
the SM prediction from the measured width, is 2 MeV (95% CL) from LEP data (52j. For h masses below the Z 
mass, we can see from Fig. [T] (a) that cos 2 (/3 — a) is always greater than 0.75 for L3. The lower bound on mt is given 
as a function of 5tan/3 in Fig. [5] for two extreme values of cos 2 (/3 — a). 

In the doublets only case (5tan/3 ~ 1), we would have had ;> 80 GeV. Since the nih and run (— 125 GeV) would 
then be of similar size, v\ and vi would be also of similar size, unless some of the X L values in the Higgs mass-squared 
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FIG. 2: Lower bound on the light Higgs mass from the Z — > hZ 1 contribution to the invisible width. Upper (Lower) curve is 
for the extreme value of cos 2 (/3 — a) = 1.0 (0.75). The doublets only (no singlet) case corresponds to 5 tan/3 ~ 1 point. 



matrix in Eq. (|35l) are larger than 4w, violating perturbativity. The requirement of perturbativity then gives a limit 
of tan j3 <J 10 (or 5 J> 0.1). As discussed in the Appendix, experimental constraints provide a tight bound on 5 in the 
doublets only case which conflicts with this limit. With a presence of the Higgs singlet, 8 tan /3 can be lower than 1, 
and the mn bound gets relaxed, as Fig. [2]shows. For Jtan/3 < 0.1, there is essentially no bound on rrih from the LEP 
invisible decay width. Thus, we find that the doublets only case is unacceptable in the scenario we consider. 

We see that for sufficiently small <5tan/3, one can have mh < 62.5 GeV, leading to the possibility of H — > hh. As 
shown in Ref . [4l[ , this decay would unacceptably dominate H decays unless sin a was infinitesimally close to 1 . We 
will not consider this possibility in this paper. 



C. LHC bounds from the 125 GeV SM-like Higgs 



The discovery of a Higgs boson, H , at 125 GeV with roughly SM branching ratios will imply stringent constraints 
on additional decays of the H. In particular, the new decays H — > ZZ' and H —> Z'Z' should be sufficiently small. 
For a sufficiently light Z' , using Eqs. ([25| and (|26|) . we have 



(6tan/3) 2 sin 2 (/3-a) (44) 



H 



and 



r( H -> Z'Z') g ^L^f Q5 tan /3) 4 f ^ ^sin a + sin 3 /3 cos a \ 
1 ' 128Trm 2 z y H ' V cos/3sin/3 / V ' 

In Fig. [31 we show the region in which their partial decay widths are within 10% of the SM total decay widths 
(4.1 McV for the 125 GeV Higgs), illustrating the sensitive parameter-space for each decay for a few values of 
#tan/3 = 0.1, 0.2, 0.3. In the region between the colored bands, the partial decay widths are more than 10% of the 
SM total decay width. A very light Z' at the LHC could be found by a narrow resonance search after appropriate cuts 
as studied in Ref. [2l],|4(J. Constraints on the model from these decay modes need a dedicated LHC data analysis. 

Figure 2] (a) shows the parameter-space where the total decay width of the heavy Higgs H is close to the SM width 
(within 20%). We have used a modified version of the code of Ref. [4l| for some of the plots in this section. It shows 
that for a relatively large 8 taxi P, the Th would be quite different from the SM prediction except for a very narrow 
region of parameter-space. 

The scenario that the heavy Higgs is the SM-like Higgs in ordinary 2HDMs were discussed in Ref. [4l|. We follow 
their parametrization to measure the event rate ratio compared to the SM predictions. For the H — > 77, 

sma\ 2 BRDg> k M (# 77) 



1 « <46) 
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-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



sin a sin a 

(a) (b) 

FIG. 3: Decays of the 125 GeV Higgs to Z and Z' , with (a) representing H ZZ' and (b) representing H ->• Z' Z' . The 
colored bands are the areas where H — > ZZ 1 , H — > Z 1 Z 1 are within 10% of the SM width. <5tan/3 = 0.1 (blue band), 0.2 (red 
band), 0.3 (green band) are illustrated. The smaller 5 tan /3 region covers the region of larger <5tan/3, i.e. the red (green) band 
is within the blue (red) band. 
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(a) (b) 

FIG. 4: (a) Region where the total Th is within 20% of the the SM prediction, (b) Region consistent with the SM predictions 
within a factor of two in all bb, t + t~ , WW, ZZ channels. Choices of 5 tan /3 and color code are the same as Fig. |3] 

with BRsM(-ff — > 77) being our calculated value, instead of the precise value of the SM [54j. Similarly, f?#(66), etc. 
are defined for other decay modes. The prefactor comes from the relative Yukawa coupling of the H and SM quarks 
[Eq. (|40p], which is relevant for the dominant gluon fusion gg — >■ H through quark loops. 

As noted in Ref. [41], the number of events (production cross section times branching ratio) of the H into many 
modes can vary up a factor of two, given currently available data. Figure [4] (b) shows the parameter-space in the 
Dark 2HDM where the branching ratios of the 66, t + t~, WW and ZZ are consistent with the SM prediction up to 
a factor of two < rju < !)■ In the region between the bands, at least one of these event rate ratios is more than a 
factor of two different from the SM prediction. 

The 77 mode is the decay mode that drove the recent discovery of the 125 GeV Higgs boson at the LHC. We plot 
the allowed region corresponding to 77^(77) > 0.8 in Fig.[5j The exact ^(77) = 1 is not achievable with given choices 
of <5tan/3 (= 0.1, 0.2, 0.3), and requires much smaller values to reach this limit. In the ordinary Type I model, it is 
well known that one can not substantially exceed f]H{ll) = 1- 

One can expand the allowed parameter-space considerably by going to smaller i5tan/3 values. We also note the 
parameter-space of the diphoton constraint in Fig. [5] is covered, except for a very little portion at the edge, by the 
parameter-space of less stringent other modes in Fig.|H(b). 
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sin a sin a 

(a) (b) 

FIG. 5: (a) Region where the diphoton rate ratio is ijh (77) > 0.8. (b) Close-up version of far-right part of the (a). Choices 
of 5 tan /3 and color code are the same as Fig. [3] 

The effect of the charged Higgs that can contribute to diphoton signals may depend on the parameters of the scalar 
potential, and it is ignored (which is valid when the charged Higgs is heavy enough). We will consider relatively light 
charged Higgs scenario later in this paper though, which can change the diphoton rate. It is worth mentioning that 
Vh{h) ~ 1 may not be a decisive constraint as additional vector-like leptons can increase the diphoton decay rate 
although they can also alter other decays [HI H(| ■ 



IV. DECAYS AND DETECTION OF LIGHT HIGGS, h 



The most interesting decay of the h would be h — > Z' Z' . (For simplicity, we consider only the nih <J m-z case to 
avoid h — > ZZ' decay.) 

For a sufficiently light Z' ', using Eq. (|!?9"|) . we have 

, , g 2 z m\ .4 ( cos 3 j3 cos a — sin 3 j3 sin a 



T(h -> Z'Z') ~ f (S tan pf 1 a . " . (47) 

v ' \2%-Km\ y ,! \ cos /3 sin /3 1 



We can compare this to the typically dominant decay channel, bb, 



3m%mh ( cos 2 



which gives 



T(h — > bb) 12ml ^ ( cos /3 sin /3 \ /cos a 



T{h^Z'Z') m 2 h (5tan/3) 4 V cos 3 /3 cos a - sin 3 /3 sin a / \sinP 



(49) 



Figure IS] shows the parameter region for mh — 60 and 90 GeV in which h — > Z'Z' dominates the light Higgs decay 
with 50%, 90% of the total h decay. The Z'Z' mode can dominate the light Higgs decay in a similar manner in which 
the Higgs decay to the weak vector bosons would dominate in the SM if the Higgs mass were sufficiently heavy. This 
originates from the enhancement from the longitudinal polarization of vector bosons when they are boosted. 

In accordance with expectation from Eq. (|4"9")l . Fig. [5] shows that the h — > Z'Z' dominates in a larger region of 
parameter-space with larger mh and larger 5tan/3. Especially, in the sin a ~ ±1 and sufficiently large tan/3 region, 
the h decay is almost entirely into Z'Z' as indicated by 90% dashed curves. 

We note that for a re case, in which the lighter one would be the SM-likc Higgs, the heavier one would dominantly 
decay into Z'Z'. For some discussions of a Higgs decaying into Z'Z' in different contexts, see Refs. |43l 15511 . 

The Z' will then decay into fermions with a partial decay width, when fermion masses are neglected, [39j 

T(Z' -> //) ~ |§4sl (g'vf +9A f ) m z > (50) 
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FIG. 6: Boundary where the light Higgs decay is dominated by Z' Z' channel (50%) for (a) m,h = 60 GeV and (b) mt = 90 GeV. 
Dotted (Dashed) curves are where the branching ratio is 50% (90%). Choices of <5tan/3 and color code are the same as Fig. [3] 



where gAf = —T 3 f and g' V f = T3/ — 2Qf (sin 2 9w — (e/ez) cos 9w sin 9w) ■ Nc = 3 (1) for quarks (leptons). The 
branching ratio of Z' into the c harg ed leptons (which we call x e and x^) depend on e/ez as well as nontrivial hadronic 
decays. It was shown in Ref. [56| that the individual lepton branching ratios for dark photons (i.e. a model with 
only kinetic mixing e for the interaction) vary, roughly, between 10% and 40% over most of the mass range. This 
branching ratio would change if there are other light hidden sector particles that Z' can decay into. Also, our dark Z 
(i.e. a model with both kinetic mixing and Z-Z' mass mixing parametrized by ez) has different branching ratios in 
general. We assume no displaced vertex and take x^ = 10% in the following sketchy analysis. 

The light Higgs might be detectable by a similar resonance search for the heavier Higgs decaying into Z' in Sec. lIII Cl 
although a feasibility study might be necessary. Very recently, ATLAS has looked for prompt "lepton-jets" at 7 TeV 
[57l ]. A lepton-jet is a final state consisting of collimated muons or electrons. In one of their analyses, with results 
consistent with the SM, they look for pairs of lepton-jets, each with two or more muons (the rest of the event is 
ignored). Although they do consider a specific Hidden Valley model as an example, their results will apply to our 
case. They only consider Z' masses of 300 and 500 MeV, and the results are relatively insensitive to the choice of 
masses. The resulting 95% CL upper bound on the cross section times branching ratio is 17 fb for mz< — 300 MeV 
and 19 fb for a mz< — 500 MeV. For our ballpark estimate, we will take the bound to be 20 fb without specifying 
the Z' mass. 

The SM cross section for Higgs production through gluon fusion is about 37 (30, 24) pb for the Higgs mass of 
80 (90, 100) GeV at 7 TeV 58]. In our case, the rate is multiplied by the factor (cos a/ sin /3) 2 of Eq. (gT|). As shown 
in Fig. [6l the h can decay into Z' Z' dominantly in most of the parameter-space of interest. If we ignore the tan /3 and 
sin a dependence and take BR(/i — >• Z' Z') m 1, the light Higgs h will appear as two muon-jets in the ATLAS analysis 
[53|. Taking rrih = 80 GeV, the cross section times branching ratio for pp — >• h —> Z'Z' — > 2 muon-jets is bounded as 

< 37pb >(if)^ * 20,b - < 51 » 

For relatively large tan/3, this gives sin a > 0.97 for a choice of x M = 0.1. As tan/3 approaches 1, this changes to 
sin a > 0.99. As the h mass varies, these numbers change slightly, but it is clear that one is forced into a fairly small 
area of parameter-space. Although these values of sin a might seem fine-tuned, one should note that it varies from 4 
to 7 in terms of tana as sin a varies from 0.97 to 0.99. 

Given the level of our rather sketchy estimate, we do not take these bounds too seriously, but it is indicative that 
a sophisticated analysis can potentially shrink the allowed region considerably (unless the effect was discovered). 



V. DECAYS AND DETECTION OF CHARGED HIGGS, H± 



We have noted that the Dark 2HDM is unusual in that there is no pseudoscalar Higgs boson (it becomes the 
longitudinal component of the Z') associated with the doublets. However, it certainly does have a charged scalar, 
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FIG. 7: The partial decay widths for charged Higgs as a function of the charged Higgs mass. The dashed lines are for the 
decay H + — > tb, with the upper (lower) line corresponding to tan/3 = 2 (10). The dotted lines are similar for the decay 
H + —¥ t + i>. The solid line is for H + — > hW + , which is either virtual or real, assuming ran = 80 GeV and cos 2 (/3 — a) = 0.75. 
The dot-dashed line is for H + —¥ HW + with the 125 GeV SM-like Higgs. The dominance of hW decay mode is persistent for 
all cos 2 (/3 — a) values allowed by the LEP over the mass range of interest: ran — 60 — 90 GeV. 



H ± , which is orthogonal to the longitudinal component of the . The phenomenology of the charged Higgs in the 
ordinary Type I model was recently discussed in Ref. fl6l ] . 

In the Type I model as well as in this Dark 2HDM, the coupling of the charged Higgs to fermions is suppressed by 
a factor of tan/3. (Thus B decay data, which pushes the charged Higgs mass above 300 GeV in the Type II model, is 
not stringent in the Type I model and in fact gives no bound for tan/3 > 2 [l^, Direct bounds on the charged 

Higgs mass from LEP H + H~ searches give m H ± > 79 GeV assuming the charged Higgs decays only into H + — > cs 
and t + v [60j. 

If a charged Higgs is light enough, one can look for it in top quark decays, but as ATLAS data shows, the lower 
bound on tan (3 in the Type I model is fairly weak, varying from 1 to 4 in the mass range of 90 GeV < m H ± < 140 GeV 
when the dominant decay mode is t + v [61|. For charged Higgs masses above about 180 GeV, its primary decay mode 
is generally given by H + — > tb, which is very difficult to detect. For some discussion of relatively light charged Higgs 
(m H ± £ 90 GeV) in ordinary 2HDMs, see Ref. (62|. 

In the Type I model, one can also look for the decays — > ipW^ where <p is either h, H or A. In the Dark 
2HDM considered in this paper, however, there is a new and dramatic signature. We have seen the primary decay 
mode of the h in most of the relevant parameter-space of our interest is into Z' Z 1 . In this case, one can have 
i7 ± — > hW^ — > Z'Z'W^. Much of the time, the Z' will decay into lepton pairs. This will certainly be the dominant 
decay if the — > hW^ decay is kinematically a two-body decay, and would still be dominant even if the decay is 
kinematically three-body, since the competing decay will be suppressed by two powers of the r mass. 

The widths for the two and three-body decays of i/ ± — > hW^ can be found in Refs. [63|,|64j]. The two-body decay 
width is 

r(i^ -> hW±) = C0S y~ Q) 3 X^ 2 (m 2 H± ,mlm 2 w ) (52) 

with A(a;, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx. It decreases with and increases with cos 2 (/3 — a). Similarly, the 
decay width into the heavy Higgs (125 GeV) is given by 



sin 2 (/3 -a) 1 3/2 2 2 2 

and the decay width into tb is given by 



r(ff± -> HW±) = ^ ~' I X 3/2 (m 2 H± , m 2 H , m 2 w ), (53) 



T(H± -+ tb) * N TT m L ( 1 - ^" ) (54) 



8ttv 2 tan 2 (3 



when rrif, is ignored and Vtb = 1 is taken, and a similar expression is given for t + v mode. 
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Figure [7] shows the partial decay widths based on the above formulas for H + — > tb, H + — > t + v, H + — > hW + , 
H + — > HW + assuming that = 80 GeV and cos 2 (/3 — a) = 0.75. (The cs has a smaller contribution than t + v as in 
the ordinary Type I model because of the small mass at the Higgs mass scale as the RG running is faster with color.) 
For uiii < m H ± < raw + rrih region, the three-body decay H + — > hW + * — > hff' (and similarly for H) is shown in 
the plot. 

We see that the decay into hW completely dominates for the entire range of masses in the plot {m H ± > 
100 GeV). We checked that, for all our interested range of mj, ~ 60 — 90 GeV and cos 2 (/3 — a) ~ 0.75 — 1 that we 
discussed in Sec. IIII1 the -ff ± — > hW^ decay mode keeps dominating over the other decay modes for the mass 
range of mjf± > 90 GeV. While the decay rate of the charged Higgs is the same as in ordinary Type I model, the 
difference comes from how the h can decay, leading to more dramatic signatures. 

The production cross section for charged Higgs bosons, for fairly large tan/3, is dominated by Drell-Yan pair 
production. Each of these charged Higgs can mainly decay into Z'Z'W ± , giving four Z's in the relevant parameter- 
space. This may be looked for with lepton-jet searches 57J in a similar fashion we discussed in the previous section 
for pp — > h — > Z' Z' channel, with appropriate cuts and selections. Note that the lepton-jet searches only require at 
least two muon-jets. With four Z's and a O(0. 1) branching ratio for Z' — > f.i + fJ.~, a substantial fraction of charged 
Higgs pairs will give a signal. One can estimate that the current ATLAS bound of 20 fb is already covering a region 
of parameter-space, but it is clear that a much more detailed analysis is needed. Given the unique signature of the 
model, with four Z's and two W's, a more targeted search could cover much more of the parameter-space. 

In the Dark 2HDM, there is no tree-level vertex H^W^Z' for the same reason that there is no vertex 
at tree level. Namely, if one goes into a basis in which only one Higgs doublet gets a vev, then the charged Higgs 
is entirely in the other doublet. As a result, it can have no vev-dependent vertices. Since SU(2)l is broken, this 
result will break down at one- loop, and thus — > W^Z and — > W Z' can occur at one- loop. The effect of 
the Goldstone boson equivalence theorem, with a boosted Z' , would enhance this loop-suppressed decay potentially 
at meaningful level. 

We discussed, in this paper, only the case of m H ± > 100 GeV. It will be interesting and important to study how 
low the charged Higgs mass can be in this model while satisfying all the experimental and theoretical constraints, in 
view of the fact that most of the lower bounds on the to#± were obtained based on the typical tb, cs and t + v modes. 
It is also noteworthy that other variants of the Dark 2HDM (for example, with some of the SM fermions carrying 
nonzero U(l)' charges) would give the similar hW dominance since Eq. (|5"2"j) is valid for all types of 2HDMs. 



VI. SUMMARY AND CONCLUSIONS 



Although experiments at the LHC have recently found the SM-like Higgs boson with a mass of 125 GeV, it is still 
important to search for other Higgs bosons, which occur in most extensions of the Standard Model. In this paper, 
we considered a 2HDM with a new U(l) gauge symmetry under which the SM particles are not charged (the Dark 
2HDM), and discussed the physics of the other Higgs bosons, namely, the nonSM-like neutral Higgs and the charged 
Higgs scalar. The additional Higgs doublet of this Dark 2HDM is charged by the U(l) gauge symmetry, and this U(l) 
plays the role of the Zi parity of the ordinary 2HDMs, thus forbidding tree-level FCNC. 

Kinetic mixing between the U(l) gauge groups will generally occur, as will mixing in the Z-Z' mass matrix. The 
possibility of a very light Z' gauge boson has attracted increasing interest of late. Such a light Z' has been a subject 
of active experimental searches including the fixed target experiments at JLab in Virginia and at Mainz in Germany. 
There are also searches using the decays from mesons at KLOE, BaBar, and Belle experiments. 

The physics connecting the Higgs to a heavy Z' (for example, the Z' decays into the Higgs boson [6^, [6(|) has been 
extensively studied, but connection of the Higgs to a very light Z' has not been. The light Z' allows a very interesting 
scenario as the various Higgs boson decays can involve a light gauge boson which can decay into leptons with 0(0.1) 
branching ratio along with possible enhancement from the longitudinal polarization. Interestingly, in this model, the 
most important predictions/constraints on Higgs properties come from the properties of the Z' . The Z' is assumed 
to have a mass of 0(1) GeV or less. Such a light mass is the region most high energy collider analyses discard to 
avoid large SM backgrounds. 

We considered the case in which the 125 GeV SM-like Higgs is the heavier scalar, and studied the light Higgs and 
charged Higgs in the Dark 2HDM. After describing the complete model, we studied constraints arising from LEP 
bounds for invisible Higgs decays and width of the Z as well as the 125 GeV SM-like Higgs state seen at the LHC. 
We found that a model with just two Higgs doublets is excluded by the LEP invisible Higgs decay bounds combined 
with the precision electroweak physics, rare B decay, etc. This issue can be resolved by introducing an additional 
Higgs singlet. 

We then considered the production and detection of the light Higgs boson (h) whose dominant decay is into Z'Z'. 
This leads to a remarkable and unusual signature. The Z' will decay with an appreciable branching fraction into a 
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lepton pair, which (since the Z' is so light) will form a collimated lepton jet. A recent ATLAS experiment looked for 
events with two or more such "muon-jets" , and their results show that significant constraints on the model would be 
possible if a sophisticated analysis with larger statistics were to follow. 

The ATLAS muon-jets experiment also constrains the charged Higgs boson, which decays predominantly into hW 
(either as a two-body or a three-body decay). The charged Higgs in the Dark 2HDM can be very light compared to 
those in the ordinary 2HDMs. 

In this paper, we limited ourselves to the case in which the 125 GeV SM-like Higgs is the heavier neutral Higgs, 
no mixing is present between the Higgs doublets and a singlet, and we focused on the region in which the light Higgs 
mass is in the range of m#/2 < w,h < mz- Relaxing these limits will allow more studies. 
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A. Hass, E. Strauss, and G. Watts for discussions about their muon-jet analysis. HL thanks H. Davoudiasl, I. Lewis, 
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Appendix A: Constraints on Z' properties 

In the m 2 z , <C m 2 z0 limit, as shown in Sec. [Hi the Z-Z' mass-squared matrix is given by 

^■ = (-hX) am% °(-(<n%Xi.)- (A1 » 

We take the parametrization £ = ez + £ tan 8w which helps separate the Z' interactions with Jnc and J em . The e 
comes from the kinetic mixing term, and the ez originates from the Z-Z 1 mass mixing from the Higgs ($1) that is 
charged under both SU(2) L x U(1)y and U(l)'. 
The mass-squared matrix then can be written as 

M ^'^ m ^{-(s z +eUn9 w ) m%,/m% ) (A2) 

taking m 2 z , ~ m z , which is realized for £ 2 <C m 2 z , /m 2 z . The Z-Z' mass mixing parameter ez is further parametrized 
by 

ez = ^6 (A3) 
m z 



(A4) 



with, from Eqs. JT3]) - (fl5j) . 

cos (3 cos (3d 



yl — cos 2 (3 cos 2 fid 



This would have been 5 ~ 1/ tan/3 in the doublets only case (cos (3d = !)■ 

The approximation 5 rs cos [3 cos (3d ~ 1/ (tan /3 tan /3d) would be valid in the limit where both tan/3 and tan (3d are 
large, and partly because of this reason, we use 5 tan /3 as our input instead of 5 in the numerical analysis of this 
paper, which helps in estimating the Higgs singlet contribution. 

Ignoring the higher order terms for small e and ez parameters (and their combination), we get 

A„t = -eJ^Ap-gzJ&cZ" (A5) 
~ -eJ? m (A, + eZ'°) ~ gzJ^ciK - £tail 0wZ'°) (A6) 
~ -eJ^Ap + EZ'J-gzJ&ciZp + ezZ'J (A7) 



where Eq. (IA6[) is obtained after field redefinition to remove kinetic mixing term [Eq. |T])] at leading order and Eq. (|A7[) 
is after the Z-Z' mass-squared matrix diagonalization [Eq. (|12p]. We can see that, because of cancellation, there is 
no net Z' coupling to the weak neutral current induced by the kinetic mixing (e). 

For a very light Z' , the kinetic mixing e is constrained by various experiments including the electron beam dump, 
electron anomalous magnetic moment, and narrow resonance searches. In particular, there are active searches using 
fixed target experiments at JLab (in Virginia) and at Mainz (in Germany) as well as the searches using the heavy 
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meson decays at KLOE, BaBar, Belle experiments. Very roughly, for the 10 MeV <; mz> < 1 GeV range, the 
experimental bounds of e 2 <; 10 -5 are present. With some combination of e and mz> (roughly, e 2 ~ 10~ 6 — 10~ 5 
and mz> ~ 20 — 50 MeV), the Z' can explain the 3.6<r deviation of the muon g — 2. (See Ref. [3!| for details.) The 
e basically parametrizes the vector coupling of the Z' as it dominates the coupling to the electromagnetic current. A 
good summar y of the constraints and sensitivities from various experiments in the e 2 — mz' parameter-space can be 
found in Ref. pl\. 

The Z' interaction to the weak neutral current is enabled by the Z-Z' mass mixing Sz, and it controls the axial 
coupling of the Z' as it dominates the coupling to the weak neutral current. The axial couplings expand the phe- 
nomenology of the light Z' into more areas including the low energy parity violation and the enhancement of the 
production for the boosted Z 1 . They were studied in Refs. [U HH, [4(| ■ They include low-energy parity violation, 
rare K decays, and rare B decays. Also, bounds from the Higgs decay exists as the presence of the H — > ZZ' and 
H — > Z'Z' modes should be still consistent with the LHC data which is consistent with the 125 GeV SM Higgs 
property. Typically, 6 < 10~ 2 — 10~ 3 , depending on various conditions, is expected to satisfy all the constraints. 
There are some caveats and dependencies on other parameters about this bounds, and we will take 8 <J 10 -2 as a 
firm upper bound that we should satisfy in this paper. This is quite a small quantity that cannot be achieved with 
the pure doublets case as l/tan/3 cannot be too small (as discussed in Sec. IIIIBj) . which is an important reason that 
this model needs a Higgs singlet. 
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